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The energies of Andreev bound states (ABS) forming in a IV-terminal junction are affected by N - 1 inde¬ 
pendent macroscopic phase differences between superconducting leads and can be regarded as energy bands in 
N - 1 periodic solid owing to the 2 n periodicity in all phases. We investigate the singularities and peculiarities 
of the resulting ABS spectrum combining phenomenological and analytical methods and illustrating with the 
numerical results. We pay special attention on spin-orbit (SO) effects. We consider Weyl singularities with a 
conical spectrum that are situated at zero energy in the absence of SO interaction. We show that the SO inter¬ 
action splits the spectrum in spin like a Zeeman field would do. The singularity is preserved while departed 
from zero energy. With SO interaction, points of zero-energy form an N - 2 dimensional manifold in N - 1 
dimensional space of phases, while this dimension is N - 3 in the absence of SO interaction. The singularities 
of other type are situated near the superconducting gap edge. In the absence (presence) of SO interaction, the 
ABS spectrum at the gap edge is mathematically analogues to that at zero energy in the presence (absence) of 
SO interaction. We demonstrate that the gap edge touching (GET) points of the spectrum in principle form 
N - 2 (N - 3) dimensional manifold when the SO interaction is absent (present). Certain symmetry lines in the 
Brillouin zone of the phases are exceptional from this rule, and GET there should be considered separately. We 
derive and study the effective Hamiltonians for all the singularities under consideration. 

PACS numbers: 74.45. +C, 85.25.Cp, 03.65.Vf, 71.70.Ej 


I. INTRODUCTION 

Superconducting junctions give rise to many interesting and 
unique physical phenomena, this being a base of the numer¬ 
ous applications in the field of quantum devices. A con¬ 
ventional Josephson junction with two superconducting leads 
hosts the Andreev bound states (ABS), that carry the super¬ 
current determined by the difference of macroscopic phases 
of the lead si. The properties of ABS may be altered by con¬ 
necting the superconductors with special materials. For ex¬ 
ample, the exchange field in a ferromagnetic junction splits 
the ABS energies in spin. This may result in the 7t-state!. 
Recent studies address topologically protected bound states 
with zero-energy, called Majorana bound states, that occur 
in the ID semiconductor nanowire with spin-orbit (SO) in¬ 
teraction, Zeeman splitting, and proximity-induced supercon¬ 
ducting gap!!. The presence of Majorana bound states in the 
junction may double the period of current-phase relation!!. 
The coexistence of SO interaction and Zeeman effect breaks 
the spin-rotation and time-reversal symmetries. With symme¬ 
try broken, the Josephson current is not an odd function of 
the phase difference. This is called the anomalous Josephson 
effect!!—. 

The Josephson junctions involving various materials have 
been mostly investigated in two-terminal setups. There is 
a recent interest in multi-terminal Josephson junctions!!—. 
Such junctions have been realized, for instance, with crossed 
InSb/As nanowires!!, where SO interaction is strong. Multi¬ 
terminal Josephson junction with N superconducting leads is 
affected by N - 1 independent phase differences. The energies 
of ABS are In periodic in all phase differences. The system 
of energy levels of ABS can be regarded as a band structure 
in a N - 1 dimensions. The phase differences play the role 
of quasimomenta!!’. The multi-terminal junctions may exhibit 


topological properties even if the superconducting leads and 
the connecting region are not made from topological or other 
exotic materials. In the case of two-terminal junction, the An¬ 
dreev levels touch zero energy only when the transmission co¬ 
efficient of the normal region is unity and the phase difference 
is ip — ±n. For the multi-terminal junctions, the Andreev lev¬ 
els can reach zero energy at some isolated points in N - 1 
dimensional space of phase differences!!!?,. Such points are 
topologically protected being the Weyl singularities studied 
theoretically in 3D solids!!. 

The energy gap closes at a Weyl point and satisfies a coni¬ 
cal dispersion in the vicinity of the point. The Berry curvature 
field is divergent at the Weyl points. They can be regarded 
as Dirac monopoles of the Berry curvature field bearing the 
topological charge, ±1. A band structure with Weyl points 
can be continuously transformed to that without the points if 
two Weyl points with opposite topological charge meet each 
other to annihilate. In a 3D solid, the SO interaction and the 
inversion-symmetry breaking are essential for occurrence sta¬ 
ble Weyl point!!!!. Very recently, the angle-resolved photoe¬ 
mission spectroscopy experiments confirmed the existence of 
the Weyl points in the 3D solids, such as TaAs? 2 24 , TaP—, and 
NbAs!!. 

Riwar et al. — demonstrate the presence of Weyl points in 
the multi-terminal short Josephson junctions. In contrast to 
solids, the Weyl point does not require SO interaction. The 
Andreev levels are Kramers-degenerate for the whole (N - 1)- 
dimensional i^-space. They discuss the transconductance to 
detect the Chern number due to the Weyl points by using one 
phase as a control parameter to switch the topological state. 
To reveal experimental signatures of the topology associated 
with Weyl points, the authors of Ref. [Hill propose the following 
scheme. They consider a 2D band structure that depends on 
the two phases, ip\ and ip 2 - The property of this band struc- 
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ture can be tuned by the remaining phase, 933 . The 2D band 
structure is characterized by a Chern number that is propor¬ 
tional to the flux of Berry curvature field through a (<^>i, ^ 2 ) 2D 
plane. The Chern number is changed by one anytime the plane 
crosses the position of Weyl singularity. The Chern number is 
observed as a quantized transconductance between the leads 
one and two, in similarity with the quantum Hall effect. 

Hech et al.~ and Padurariu et air 1 also study the Andreev 
levels at energies close to zero in three-terminal Josephson 
junction. The authors claim that the zero energy states in such 
junctions may open opportunity for a single fermion manip¬ 
ulation. In the three-terminal junctions, the Weyl singularity 
is generally absent although the energy of the ABS can pass 
zero. The authors obtain a condition for zero energy ABS and 
study the density of states in detail. When SO interaction is 
present, the energy levels of the ABS are split in spin. 

In this study, we investigate theoretically the singularities 
of the ABS spectrum in four-terminal Josephson junction tak¬ 
ing SO interaction into account. We thus attempt to formulate 
the full picture of such singularities combining phenomeno¬ 
logical and analytical methods and illustrating it with numeri¬ 
cal results. The ABS energies are found from the Beenakker’s 
determinant equation^ using the scattering matrix of the junc¬ 
tion. Mostly we concentrate on the case of short junction 
where we can disregard the energy dependence of scattering 
matrix. Sometimes the absence of energy dependence leads 
to extra degeneracy in the spectrum. To lift those we take into 
account the energy dependence by perturbation theory. 

Firstly, we concentrate on Weyl singularities that occur at 
£ = 0 in the absence of SO interaction. We consider the 
behavior of the singularities upon gradual increase of the 
strength of SO interaction. We have found that SO interaction 
splits the spectrum of ABS in spin. Very much like Zeeman 
effect would do. The conical points are departed from E = 0 
to mirror symmetric positive and negative energies. The Weyl 
singularities thus remain topologically protected. A small 
modification of scattering matrix by a parameter changes the 
position of the mirror symmetric conical points, not eliminat¬ 
ing them. As we show in numerical illustration, an annihila¬ 
tion of Weyl points of opposite charge can take place upon the 
tuning of the scattering matrix. We derive the effective Hamil¬ 
tonian describing the vicinity of the Weyl points. The cones 
in the vicinity intersect E = 0 at a 2D surface in a 3D space of 
the phases. We prove that this is the general property of ABS 
spectrum in the presence of SO interaction. 

Singularities of a different type arise when ABS energies 
approaches the gap edge, E = A. Owing to mirror symme¬ 
try of Andreev spectrum, there is a level with E = -A at 
the same position in the space of the phases. We formulate 
a mathematical analogy that permits to map Weyl singulari¬ 
ties at £ = 0 in the presence (absence) of the SO interaction 
to singularities at £ = A in the absence (presence) of the SO 
interaction. Since the ABS energies in the presence of SO in¬ 
teraction reach zero energy at 2D surface, we expect the gap 
touching point to form 2D surfaces in the absence of SO inter¬ 
action. Indeed this can be seen in a concrete numerical calcu¬ 
lation. Employing the same analogy from the fact that without 
SO interaction the ABS energies reach zero at isolated point 


only, we derive that the gap edge touching (GET) in the pres¬ 
ence of SO interaction generally occurs only at the isolated 
point. This implies that even a weak SO interaction removes 
the GET. We construct the effective Hamiltonian to describe 
this situation. In this case, a weak energy dependence of scat¬ 
tering matrix can also become important. An important pecu¬ 
liarity of the ABS spectrum concerns the vicinity of symmetry 
lines in 3D elementary cell of the space of the phases. Three 
of the four superconducting phases are the same at a symmetry 
line. Therefore the four-terminal junction at a symmetry line 
can be regarded as a two-terminal junction with unequal num¬ 
ber of conduction channels in the two leads. As mentioned in 
Ref. |28| for two-terminal short junction, the SO interaction is 
irrelevant to causing spin splitting. Thus the vicinity of sym¬ 
metry lines requires a separate consideration. We derive an 
effective Hamiltonian to incorporate the details of the GET in 
the vicinity of the symmetry lines. 

These two types of singular point would reveal the topolog¬ 
ical nature of multi-terminal Josephson junctions. It brings a 
goal to propose nanostructures as artificial exotic materials. 

The structure of this article is as follows. In Sec. II, we 
explain the model for the multi-terminal Josephson junction 
and the equation to determine the ABS energies. In Sec. Ill, 
we describe the spin splitting of Weyl singularities. Section 
IV is devoted to the GET point in general. Here we formulate 
and employ the mapping between £ = 0 and £ = A. We also 
concentrate separately at the vicinity of the symmetry lines in 
this Section. We conclude in the last Section. 


II. MODEL AND FORMULATION 

In this Section, we explain the model in use. 

We consider a junction connected to N superconducting 
leads. An example of a physical system of this sort is given 
in a Fig. Q3 a )- All microscopic detail of the junction can be 
incorporated into the scattering matrix l e h of the electrons 
and holes (Fig.Q]b)). We assume certain numbers A, of spin- 
degenerate transport channels in each lead, so s e ,h are 2Mx2M 
matrices, M = Emi* A/. The symmetry of Bogoliubov de- 
Gennes equation implies the relation ,v h (£) = -gs*(-E)g with 
g = -i&y being a matrix realizing a time-inversion in the spin 
space. In addition to this, the requirement of time reversibil¬ 
ity implies s e (E) = -gsJ(E)g. The superconducting leads do 
not provide extra potential scattering. They are described by 
the Andreev reflection amplitudes for converting electron to 
hole and hole to electron that do not change transport channel 
index and are presented by diagonal matrix rhe.eh = 
ipk being a superconducting phase of a lead corresponding to 
the channel k. Xk is an energy-dependent phase of Andreev 
reflection. We assume the same material for all of the super¬ 
conducting leads, so that their order parameters are the same 
A, = A. In this case, the energy-dependent phase is the same 
for the all leads and is given by^ = arccos(£/A). The eigen¬ 
vectors of the ABS satisfy ij/ = /'chAhUie V<A- Therefore the 
energies of ABS are determined from=2 

det (e i2x — S(0, £)) = 0 


(1) 
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FIG. 1: (a) Schematics of a four-terminal junction based on crossing 
semiconductor nanowires, (nanocross)^. (b) Scattering model for 
multi-terminal Josephson junction. The junction itself is described 
by the scattering matrix S. The superconducting leads provide An¬ 
dreev reflection with amplitude Iar. We assume the same order pa¬ 
rameter in all superconducting leads while the phases can 

differ. By virtue of gauge invariance, one of the phases can be con¬ 
veniently set to zero. 

with 

S(0,E) = -gs*(0,-E)gs(0,E), (2) 

s(0,E) = e- i<t>l2 UE)e +i * 12 . (3) 

The hat symbol on s and S is omitted for simplification. Note 
that S (0, E) is a unitary matrix. 

For numerical calculations, the scattering matrix s e is taken 
as a random matrix. The SO interaction is taken into account 
in s e . In the presence (absence) of SO interaction, random 
matrix is a member of symplectic (orthogonal) ensemble. We 
introduce a parameter pso tuning the strength of SO interac¬ 
tion^. The parameter varies from 0 to 1, providing a con¬ 
tinuous transition between the orthogonal (pso = 0 ) and the 
symplectic (pso = 1 ) ensembles. 

Let us discuss the ABS energies in the vicinities of zero 
energy, E = 0, and superconducting gap edge, E = A, where 
the first term in the determinant in Eq. (QJ becomes e' 2x - -1 
and +1, respectively. In the absence of SO interaction, the 
scattering matrix s(0, E) commutes with g. Thus, 

5 (0, E) = s*(0, -E)s(0, E). (4) 

When the SO interaction is present, we can introduce a unitary 
matrix u{0, E) = gs(0, E). Since g is a real matrix, 

S (0, E) = -u(0, —E)u(0, E). (5) 

Let us disregard energy dependence of scattering matrices. 
Suppose that for a sufficiently general unitary matrix v, there 
is eigenvector \i//) = Such eigenvector would corre¬ 

spond to an ABS at E - 0 in the presence of SO interaction 
and to an ABS at E = A in the absence of SO interaction. 


Correspondingly, an eigenvector -|t f/) = v*v\4>) would be an 
ABS at £ = 0 in the absence of SO interaction and to an ABS 
at E — A in the presence of SO interaction. We thus establish 
a mapping between the situation at E — A (E = 0) with SO 
interaction and at E = 0 (E — A) without SO interaction. 


in. WEYL SINGULARITIES AT E * 0 

In this Section, we concentrate on the Weyl singularities. In 
our numerical illustrations, we concentrate on four-terminal 
short junctions with a single channel in each lead, N, = 1 and 
M — 4. In this case, the ABS energies are periodic function 
of three independent phases 0 = so we can re¬ 

strict ourselves to the Brillouin zone \g>j\ < n for i = 1,2,3. 
The time-reversibility manifests itself as the inversion sym¬ 
metry in this Brillouin zone, E(0) = E(—0). To start with, we 
demonstrate that in the absence of SO interaction the energy 
levels of ABS exhibit the band gap closing points at zero en¬ 
ergy^. Next, we continuously change the scattering matrix 
increasing the parameter pso, that is strength of SO interac¬ 
tion. We demonstrate that the SO interaction splits the conical 
spectrum in spin. The conical point departs from E - 0 and 
the energy levels cross zero energy at 2D surface rather than 
isolated point. This proves topological protection of the Weyl 
singularity. To prove the generality of our conclusions, we de¬ 
rive an effective Hamiltonian from Eq. (0} that is valid in the 
vicinity of the singularity and at weak SO interaction. 


A. Energies of the Andreev bound state 

We obtain the ABS energies from Eq. (jT|). The scattering 
matrix s e is random. In the absence of SO interaction, we 
chose the random matrix from the circular orthogonal ensem¬ 
ble and disregard its energy-dependence. 

We examine the spectrum for many random matrices. 
About 6% of them show a gap closing of the ABS energies in¬ 
dicating the Weyl singularities. The singularities always come 
in groups of four. Figure |2a) gives the positions of singular¬ 
ities for a random matrix of choice. The time-reversal invari¬ 
ance guarantees that the singularity at 0 ^ = (sc , i° > , <pf\ <pf *) is 
accompanied by the singularity of the same topological charge 
at -0 (O) . Owing to this, the positions of four singularities and 
the center of Brillouin zone (0,0,0) are in the same 2D plane. 
Figure [2lb) shows the lowest positive energy of the ABS in 
this plane. The spectrum is symmetric with respect to phase 
inversion 0 —> -0. In the origin, this energy reaches maxi¬ 
mum E = A. It drops down to two valleys close to the edges 
of the Brillouin zone. A zoom into a valley (Fig.[2jc)) shows 
that the energy actually reaches zero in two isolated points. 

In Figs. [2) d) and (e), we show the ABS energies versus ip\ 
at fixed ip 2 , 3 . In Fig.[2jd), the choice of <^> 2,3 is such that a sin¬ 
gularity is reached at some tp\. In Fig. [2}ej, the line spanned 
by changing p\ passes close to a singularity. We see the spec¬ 
trum is conical near the singularity. For A, = 1 in all four 
leads, they are four positive and four negative ABS energies. 
Since the SO interaction is absent, the energies of ABS are 
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FIG. 2: (Color online) Energies of the ABS in the four-terminal 
junction for a particular choice of S e . The SO interaction is absent 
(Pso = 0). For this choice, there are four Weyl singularities at E = 0. 
(a) The positions of the Weyl singularities in the 3D space of the 
phases. All four points and the origin lie in the same 2D plane, (b) 
The gray scale plot of the lowest positive ABS energy in the 2D plane 
shown in (a). Contour lines indicate E/A = 0.2, 0.4, 0.6, and 0.8. (c) 
A zoom of the energy landscape (b) in the vicinity of Weyl singulari¬ 
ties. The dots give the singularity positions. The dashed contour line 
indicates E/A = 0.025. The solid contour lines corresponds to E/A 
in multiples of 0.05. Panels (d) and (e) show the ABS energies when 
two of the three phases are fixed to (d) ip 2 — 0.53;r and ip 2 - -0.42 n 
and (e) ip 2 — 0.53tr and ip 3 = -0.607T. 


doubly degenerate. In Figs.[2jd) and (e), the second ABS band 
is close to +A. 

At the same choice of random scattering matrix, we in¬ 
crease the parameter pso thereby together continuously in¬ 
creasing the strength of SO interaction. As we see in Fig. 
[3a), the SO interaction splits the ABS energies in spin. The 
absence of time-reversibility required for such splitting comes 
about non zero ( 01 , 2 , 3 . As in Fig. [2] we plot in Figs. [3b) and (c) 
the lowest positive energy of the Andreev states in the plane 
that passes through the Weyl singularities. We see in Fig. [3b) 
that overall energy landscape have not changed significantly 
in comparison with Fig.[3b). However, as seen in Fig. [3c), 
the landscape has changed drastically in the vicinity of the sin¬ 
gularities. The gap is closed at the closed contour encircling 
the singularities. The singularities are shifted to non zero en¬ 
ergy while the spectrum remains conical in the vicinity of a 
singularity (Fig. [3a)). In Fig. [3d), we plot energy difference 
between the second and the lowest positive energy levels and 
observed that it goes to zero at the position of the singularities. 



FIG. 3: (Color online) Energies of the ABS in the same junction with 
SO interaction (p so = 0.3). (a) The ABS energies at ip 2 — 0.62 n and 
tpi — —0.38zr versus ipi (c. f. Fig. [3d)). All the levels are split in spin. 
The levels pass the position of the singularity at this choice. The 
singularity is shifted from zero energy while the spectmm remains 
conical at this point. The levels cross zero energy at both side of 
singularity, (b) The gray scale plot of the lowest positive ABS energy 
in the 2D plane that includes the singularities. Contour lines are the 
same in Fig. [3b), and the overall energy landscape does not change 
significantly. The qualitative change in the zoom view of (c). Thick 
solid line indicates zero energy. The dots indicate the singularities 
which are now local maximums of the lowest positive ABS energy. 
Thin solid contour lines give E/A in multiples of 0.05. The dashed 
line indicates E/A = 0.075. (d) The gray scale plot of the energy 
difference between the second and the lowest positive ABS energies. 
The contour lines indicate E/A in multiples of 0.05. The energy 
difference reaches zero at the positions of singularities. 


The SO interaction changes the position of the Weyl sin¬ 
gularities while preserving their topological charge and the 
conical dispersion. When pso increases gradually, the four 
Weyl singularities move in the 3D space of the phases. Fig¬ 
ure [3 a) shows trajectories of their positions for pso arrang¬ 
ing from pso - 0 to pso — 0.462. Solid and dashed curves 
give the position of the singularities with positive and nega¬ 
tive charge. For this particular choice of the scattering matrix, 
the singularities of the opposite charge get close to each other 
upon increasing the SO strength and eventually annihilate at 
Pso — 0.462, so the junction is not topological any more (Fig. 
0a)). 

We compute the energy of the singularity and the plots re¬ 
sult in Fig. 0b). We see that this energy is zero in the absence 
of SO interaction. The energies of the singularities of differ¬ 
ent topological charge increase and become different with in¬ 
creasing pso- At pso > 0.3, the energies come close together 
merging at the annihilation point pso — 0,462. 


B. Effective Hamiltonian 

To prove the generality of the above results, we derive an 
effective Hamiltonian for ABS that is valid in the vicinity of a 
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For 3D space of the phases, the matrix M is positively de¬ 
fined. We reproduce a conical spectrum of ABS in the vicin¬ 
ity of singularity. For N dimensional space, the matrix M has 
N - 3 zero eigenvalues. We stay at zero energy if we depart 
from ( f (0) in this N - 3 directions. So that, the singularities are 
concentrated at N - 3 dimensional manifold. 

Let us take into account weak SO interaction, expanding 


FIG. 4: (Color online) (a) The trajectories of the positions of Weyl 
singularities in a 3D space of the phases upon increasing the SO 
strength from pso = 0 to pso — 0.462 when the singularities of 
the opposite topological charge come together and annihilate. Solid 
(dashed) curves give the position of the singularities with positive 
(negative) charge, (b) Energies of the singularities with positive 
(solid) negative (dashed) topological charge versus p so . The ener¬ 
gies merge in the annihilation points. 


singularity. When the SO interaction is absent, the singularity 
is at £ = 0. At the position of singularity, 0 (O \ the scatter¬ 
ing matrix S (0 {)) , 0) should have an eigenvalue -1. However, 
S is a rather special matrix: at £ = 0 it can be presented as 
S = s ( 0 ) *s (0) (cf. Eq. ©). Here we introduce s (0) = s(ip ( 0 ) ,0). 
We assume nothing about ,v (0) regarding it as an arbitrary uni¬ 
tary matrix. This implies that the eigenvalues of S come in 
complex-conjugated pairs. If a) is an eigenvector of S with 
eigenvalue e a , \b) = i ( 0 , *|a*) is an eigenvector with eigen¬ 
value e~ a . This implies that the two orthogonal eigenvectors 
corresponding to eigenvalue - 1 , 

s (0) *5 (( Vl?> = -\o,b). (6) 


To obtain the effective Hamiltonian, we project the matrix in 
Eq. O on to the space spanned by these two eigenvectors. For 
small deviations of the phases from < p (0 \ 50 - ip - 0/ O) <s 1 . 
We expand the scattering matrix as s( 0) = s^e' x( - s ® ~ s (0) (l + 
iX), where X is a Hermitian matrix proportional to 60. Up to 
the first order in X and £, we find 

S(0) * s (0) V 0) + iX; X= s (0) *s (0> X-s (0> X*s (0 \ (7) 
e i2x * -1 + t'2£/A. ( 8 ) 


With this, the determinant equation (jT|) can be presented as an 
eigenvalue equation for a 2 x 2 effective Hamiltonian^, 


E = H; 


H = A £*;£/ = 
7=1 


t-i 11 2 


{a\X\a) {a\X\b) \ 
(b\X\a) (b\X\b) ) 


(9) 


where £ 1 , 2,3 are the Pauli matrices in the basis of |n) and | b). 
Using the relation | b*) = / 0 ) |a), we prove that (b\X\b) = 
~(a\X\a) and real parameters ^ 1 , 2,3 are given simply by X\ + 
iX 2 = -2{b\X\a) and A 3 = -(a\X\a) + (b\X\b). The Hamilto¬ 
nian in Eq. © is a case of Weyl Hamiltonian. 

It has opposite eigenvalues, £ = +A ^X 2 + X 2 + X 2 . Ex¬ 
panding in 5p, Xj = X' m 5(p m , we obtain 



s(^ 0) ) - „ s (0) 

S {0 O) ) = -gs*(^ 0 ) )gs(^ 0) ) » s 



( 0 )* 40 ) 






( 11 ) 


K a = s' 


,( 0 )« ( 0 ) 


K a + s^KU 


dO) 


( 12 ) 


with fr a being Pauli matrices in spin space, K„ being associ¬ 
ated Hermitian matrices in channel space. We project on four 
dimensional space of spins and vectors | a, b). We observe that 
the structure of matrices K a is quit difference from X. Using 
the relations between | a) and | b), we prove (a\K a \b) = 0 and 
(a\K a \a) - (b\K a \b) = ~{a\K a \a) + (b\K a \b) = Kq. With this, 
the effective Hamiltonian becomes: 

3 

H = A ^ x j^j + A& ' K o- (13) 

7=1 

The BdG symmetry thus guarantees a special structure of this 
Hamiltonian where spin and orbital degree of freedom are to¬ 
tally separated. Kq plays a role of an effective Zeeman field 
that splits the original conical spectrum. 


E/A = cr k \Kq\ ± 



(14) 


with ctk = ± being the spin projection on the axis of the ef¬ 
fective Zeeman field. The SO interaction does not remove 
the conical point but rather shift it in energy by ±A|£'o|. In a 
3D space, the zero energy is achieved at 2D surface of ellip¬ 
soid defined by Y.k,m ^VmMmkbpk = | A'o1 2 - The singularity is 
enclosed by the surface. This consideration shows generality 
of our numerical results. The size of the ellipsoid enclosing 
a singularity increases with increasing SO interaction. In Fig. 
[3jc), we see that the ellipsoids enclosing each singularity have 
already merged together at pso = 0.3. 


IV. GAP EDGE TOUCHING 


In this Section, we consider the ABS in the vicinity of the 
superconducting gap edge £ « A. We show that the ABS 
energies reach the GET in the absence of SO interaction at a 
2D surface in the 3D space of the phases. The SO interaction 
lifts the GET almost everywhere except particular manifolds: 
symmetry lines and isolated points. We investigate these cases 
separately and establish effective Hamiltonians. 
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A. GET at symmetry lines 

To understand the peculiarities of the GET for multi¬ 
terminal junction, let us first consider two-terminal one with 
unequal number channels in the left and right lead, Nl < Nr. 
The estimation of a number of localized Andreev states is 
somehow ambiguous. From one hand, one may argue that 
there are only Nl such state because only that many states are 
sensitive to the superconducting phase difference between two 
leads. From the other hand, the full number of Andreev states 
is given by (Nl + Nr)/2 provided the total number of chan¬ 
nels is even. The remaining (Nr - AjJ/2 states are pinned to 
the gap edge with no regard for superconducting phase dif¬ 
ference. In a two-terminal junction, these states are indistin¬ 
guishable from the states of the continuous spectrum. This 
is not a case of multi-terminal junction. We note that a multi¬ 
terminal junction in fact becomes two-terminal one if the leads 
are separated in two groups with superconducting phases are 
the same within each group. For instance, in our four-terminal 
set up, one can choose <po — 0 and ipi — (p 2 — <Pi + 0. Such 
setting defines a symmetry line in multi-dimensional space of 
the phases. For our example with one channel in each lead, 
we find extra state pinned at the gap edge along the symmetry 
line. In distinction for a two-terminal case, this extra state can 
not be attributed to the continuous spectrum since it departs 
from the gap edge if we go off the symmetry line. 

Since the SO interaction does not work for a short two- 
terminal junction^, the SO splitting is also absent at these 
lines. We postpone the discussion of the details of the spec¬ 
trum in the vicinity of symmetry lines to Sec. II V Dl and con¬ 
centrate now on general situation. 

B. ABS energies near the gap edge: general 

Fet us consider GET at general position in 3D space of 
phases. Fet us note mathematical analogy between the spec¬ 
trum at E — 0 and E — A. The spectrum is determined by 
properties of the scattering matrix S in Eq. ©. Zero energies 
(E = 0) correspond to the eigenvalue of -1 of the matrix while 
GET correspond to the eigenvalue of 1. From the other hand, 
as we have seen in Section HU in the absence and presence of 
SO interaction, we can represent as S = v*v and S = —v*v, re¬ 
spectively, v» being a general unitary matrix. This establishes 
a rather unexpected mapping of the spectrum at £ = 0 in the 
presence (absence) of the SO interaction to the spectrum at 
E = A in the absence (presence) of the SO interaction. We 
have derived in Section ITlI Al that the ABS reach E = 0 at 0D 
or 2D manifolds in the 3D spaces in the absence or presence 
of SO interaction, respectively. This implies that the GET oc¬ 
curs at 2D and 0D manifolds in the absence and presence of 
SO interaction, respectively. 

Fet us see this in numerical results. In the absence of SO 
interaction, scattering matrix s e is chosen randomly from the 
circular orthogonal ensemble. Its energy-dependence is disre¬ 
garded. We examine the spectrum for many random scatter¬ 
ing matrices. The GET is observed for all the matrices. For 
a 3D space of the phases, the ABS energy reaches E = A 



FIG. 5: (Color online) GET in a four-terminal junction in the absence 
of SO interaction, (a) Curve of the GET at fixed = -0.8rr in the 
3D space of the phases. Four solid lines are symmetry lines, where 
three of the four superconducting phases are equal, (b) Curves of the 
GET at i />3 = -n (left upper), —0.8 tt (left middle), -0.6 n (left bot¬ 
tom), -0.4 tt (right upper), -0.2 n (right middle), and 0 (right bottom 
panels). Dashed curves in right bottom panel indicate if 3 = -0. In. 
The black dots give intersection points of the symmetry lines with 
the corresponding ipi-ifi 2 planes. 


at a 2D surface. This is illustrated in Fig. [5] We fix tf \ and 
find the GET points forming a ID curve. Interestingly, the 
curve passes symmetry lines where we indeed expect GET. 
For instance, in Fig. 0 a), the curve passes the symmetry line 
(pi — <p 2 = (f 3 at ((pi, (pi) ~ (— 7 T, — n) and passes the symmetry 
line (pi = (pi — 0(= tpo) at ((pi,(pi) = (0,0). As <^3 is tuned 
from - 7 r to 0 (Fig. 0b)), the positions of the first passing is 
gradually sifted to ((pi,(pi) ~ (0,0). At ^3 = 0, the curve 
approaches two symmetry lines in this plane, ipi = - 0 

and ipi = >P 3 — 0. For positive ip^, the spectrum is obtained 
from the inversion symmetry. We see that the GET occurs at 
a single connected 2D surface that includes all four symmetry 
lines. In principle, there is nothing to for bit more sophisti¬ 
cated topology of the surface. However, in a dozen of samples 
we have explored, we have found no complex topology. 

Figures0a) and (a’) show the ABS energies with the GET 
in the absence of SO interaction. The energies are plotted 
versus <pi while <p 2 and <^3 are fixed. The GET is found at 
(pi — 0.391 n. In the vicinity of the touching point, the spec¬ 
trum of ABS energy is parabolic rather than conical. A mirror 
symmetry of the ABS energies guarantees the GET at E = -A 
at the same position in the 3D space of the phases.The SO in¬ 
teraction (Figs. (b)[ 6 jind (b’)) lifts a degeneracy of the energies 
in spin. In Fig. 0b’), one spin-resolved ABS energy comes 
very close to the gap edge. However, as seen in the inset, the 
touching does not take place in the presence of SO interaction 
conform to our expectation. 

In Fig. 0 we illustrate the situation in the vicinity of sym¬ 
metry lines. We plot the ABS energies along the line <p\ = 
7t/2. This line crosses the symmetry lines. The crossing points 
are indicated by arrows in the Figure. We see that at these 
points, (p 2 — (pi = 0,7t/2, the GET survives in the presence of 
SO interaction. Also, the spin splitting of the lowest ABS van¬ 
ishes at this point. Then, this proves that the SO interaction is 
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FIG. 6: (Color online) Energies of the ABS in a four-terminal junc¬ 
tion. The sample is the same as that in Fig. [5] (a) ABS energies 
versus ip\ when spin-orbit interaction is absent, ipi and </>3 are fixed 
at 0.3/r and 0.57T, respectively. The GET is obtained at ip\ — 0.391;r. 
Panel (a’) and its inset are enlarged views, (b) ABS energies when 
spin-orbit interaction is present (p so = 0.5). The other parameters 
and scales of axes are the same as those in upper panels. 

irrelevant at the symmetry lines. 

C. effective Hamiltonian in general case 

To prove the generality our numerical results concerning 
the GET, we derive an effective Hamiltonian for ABS near 
E - A taking into account weak SO interaction and energy- 
dependence of the scattering matrix. We assume no vicinity 
of symmetry lines. This specific situation will be considered 
separately in the next Section. 

Let us first neglect SO interaction and energy-dependence 
of the scattering matrix. For a GET point at 0 ()) , the matrix 
S has to have an eigenvalue 1. We put s^ 0 ') = s (0) . This 
eigenvalue is double degenerate; if \a) is an eigenvector of 
S( v? 0 - 1 ) = s^*s^ 0) with the eigenvalue 1, | b) = s ( 0 ) |a*) is also 
an eigenvector of S (^ 0) ), 

s ( 0 ) *s ( 0 ) |a, b) = \a, b). (15) 

We project the matrix in Eq. (Q} on to together the space 
spanned by these two eigenvectors. Next we expand full scat¬ 
tering matrix in Eq. © with respect to small phase deviation 
from weak SO interaction, and weak energy-dependence. 
The first order expansion in X (deviation) and K a (SO interac¬ 
tion) is the same as given by in Eqs. (|7Ji and ( IT 2 l) . respectively. 
The energy-dependence is expanded with a Hermitian matrix 
d, 

s (0 \E ) = s <f» e iEd ^ 0)) » s (0, (l + iEd) (16) 

S(0°\ A) = s ( 0 ) *(-A)s ( 0 '(A) * s (0) *s (0) + iD; 

D = A(s ( 0 ) V°V + v (0) *cf s (0) ). (17) 

The causality of scattering matrix s(E) implies that all eigen¬ 
value of d are positive. 

The first order correction terms are similar to those for 
Weyl singularity at E » 0. However, since the sign of 




FIG. 7: (Color online) Energies of the ABS in the vicinity of symme¬ 
try lines. The sample is the same as that in Fig. [5] (a) ABS energies 
on a line tp 2 = <p 3 at ifii = n/2. A strength of spin-orbit interaction is 
pso = 0.5. (b) Enlarged view of panel (a). Arrows indicate points at 
the symmetry lines. 


the eigenvalue is opposite, the selection rules for matrix ele¬ 
ments representing the deviation and SO interaction are inter¬ 
changed. For small phase deviations, {a\X\b) = 0; (a|A|c/) = 
-(b\X\b) = (a\X\a) - (b\X\b) = A 3 . The expansion of the 
weak SO interaction gives three independent vectors in spin 
space, K 0 = « a\K\a ) + (b\K\b)) = < a\K\a) = {b\K\b > and 
K\ + 1 K 2 = 2{b\K\a) = (b\K\a). For weak energy-dependence, 
(a\D\a) = (b\D\b ) = (ci\D\a) + (b\D\b) = Z)q. No restric¬ 
tion applies to the off-diagonal matrix elements. So we define 
Z>i + 1 D 2 = 2(b\D\a) = ( b\D\a ). The positivity of d implies 

A)> ^D\ + D\. 

To deviate an effective Hamiltonian, we introduce the en¬ 
ergy deviation from the edge e — 1 - \E\/A. Up to the first 
order in -\fe, we obtain e l2x » 1 ± i V8e. Summarizing all 
these, we obtain that the spectrum is defined by the following 
eigenvalue equation 

2 2 

H = V2e; H = X 3 i 3 + ^ d- ■ Kjtj + ^ Djlj. (18) 

2=0 2=0 

Here £ 0 = 1. The full analysis of this Hamiltonian is involved. 
Here, we discuss several simple cases. Let us first neglect 
both SO interaction and energy-dependence. The eigenvalues 
of the Hamiltonian are in this case simply +IA 3 I. The nega¬ 
tive eigenvalue does not lead to any localized state. The posi¬ 
tive eigenvalue provides e = Xj/2 since A 3 is linear in phase 
deviation this defines a parabolic spectrum touching the gap 
edge at A 3 = 0. We see that the GET requires one condi¬ 
tion A 3 = 0 to be fulfilled in distinction from three condition 
A j = 0 required for Weyl singularity. This is why the manifold 
of GET point in A-dimensional space of phases generally has 
dimension N - 1 in the absence of SO interaction and energy- 
dependence of scattering matrix. 

If we take energy-dependence into account, the eigenvalues 

are given by Z>o ± X 2 + D 2 + /);. The plot of eigenvalues 
and spectrum is given in Fig. [ 8 ja). We see that the energy- 
dependence modifies the GET in a rather complex way. The 
spectrum in the vicinity of the GET line forms two bands. 
While one of the band never touches the edge, the other band 
exist only in the vicinity of the line and merges with the con¬ 
tinue upon increasing A3. 
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FIG. 8: (Color online) Fine structure of GET. X 3 gives the distance 
from the GET. (a) The effect of energy-dependence. We set D 0 = 0.3 

and yjD] + D\ =0.1. (b) The effect of SO interaction. K 0 is parallel 

to e 3 . We set K 0 = 0.1 and y/\Ki \ 2 + l/Gp = 0.2. Directions of the 
vectors are chosen randomly, (c) The combined effect of the energy- 
dependence and SO interaction. Dj and Kj are the same as those in 
(a) and (b), respectively, (d) Weyl singularity near the gap edge. We 
put Xi = A 2 = 0 and D 0 = 0.3. 


D\ {a\D\b) = 0. The effective Hamiltonian thus reads 

3 

H = Yj X £j + A>£o. (21) 

j= 1 

with X\ + iX 2 = (b\X\a) = 2{b\X\a), Xj = {a\X\a) = 
~{b\X\b) = (a\X\a) - (b\X\b), and O 0 = (a\D\a) = (b\D\b ) = 
{a\D\a) + {b\D\b) > 0. We keep here energy-dependent term 
oc Dq that is absent in the Hamiltonian ([o]i. We have just found 
a Weyl singularity with a conical point near the gap edge. The 
conical point requires three conditions to be fulfilled X t - 0. 
Therefore we expect the points to form N - 3 dimensional 
manifold in the N dimensional spaces. The energy-dependent 
term shifts the energy of the conical point from the gap edge 
similar to the effect of the SO interaction in Eq. (IT3l > (Fig. 
[8}d)). We have tried to find this singularities in our numerical 
simulation. So far, we have found none. The reason of this is 
not completely clear for us. We hypothesize that probability 
to find a random scattering matrix with such singularity is low 
because the GET point in the presence of SO interaction tend 
to stick to the symmetry line. 


Let us neglect now the energy-dependence and take into ac¬ 
count SO interaction. When K , terms are taken into account, 
the four eigenvalues of H are 


± yjx 2 +L± 2 yf(K 0 -K\) 2 + (K 0 ■ K 2 ) 2 + \X 3 K 0 - K x x K 2 \ 2 

(19) 

with L = \Ko \ 2 + \K \\ 2 + |1G| 2 - Two positive eigenvalues de¬ 
fine two spin-split bands that never touch the edge (Fig. [ 8 jb)) 
unless a special condition discussed below is fulfilled. The 
minimal energy distance to the gap edge is not achieved at 
X 3 = 0 but rather is shifted depending on the parameters of 
the SO interaction. 

If the energy-dependent terms and the spin-splitting are the 
same order of magnitude, the picture of the GET can become 
complex. For instance for an example given in Fig. [He), there 
are only three bands. One existing only in the vicinity of the 
gap edge line. 

Let us now discuss a special condition of the GET in the 
presence of sufficiently strong SO interaction. To analyze this, 
let us derive an effective Hamiltonian with strong SO interac¬ 
tion assuming a GET point to be present at 0 I]> . Strong SO 
interaction guarantees that only two spin-dependent eigenvec¬ 
tors are important instead of the four as in the previous con¬ 
sideration. The two eigenvectors | a) and | b) = M ( 0 , [a*) satisfy 


D. The vicinity of a symmetry line 

Let us consider an effective Hamiltonian in the vicinity of 
symmetry lines. The SO interaction is assumed to be strong. 
We concentrate on a four-terminal junction. At the symmetry 
lines, three of the four superconducting phases are equal. The 
electron scattering matrix can be presented in a block struc¬ 
ture. 


n 

^13 

h 1 





Assuming there are A 3 channels in three leads having the same 
phase and N\ channels in the other lead, r 3 is a 2A 3 x 2 A 3 re¬ 
flection matrix for three leads having the same phase. r\ is 
that of the other lead, tn and 131 are the transmission matrices 
between the three leads and the other lead, their dimensions 
are 2Ai x 2 A 3 , 2 A 3 x 2Ai, correspondingly. At the symme¬ 
try line, we have 2 (N 3 — Ai) independent vectors satisfying 
follA) = 0. These states are disconnected from the other lead. 
Their energies are precisely at the gap edge, therefore 


~gs 


( 0 )* &t ( 0 ) 

S ' 5 


' 0 ' 


' 0 ' 

, l<A>, 


, 1^), 


(23) 


-u ( 0 ) *u ( 0 ) \a,b) = \a,b) (20) 

withM (0) = m(<^ 0) ). The problem is thus mathematically equiv¬ 
alent to our consideration of Weyl singularity in Sec. IHIBI 
The first order expansion in X (deviation) and A d (energy- 
dependence) results in correction terms X = —u ( 0 > *u (0) X + 
u (0)* x *u (0) and D = -A (u ( 0 ) *u <0) d + M (0) *cf m ( 0) ), respectively. 
We prove the selection rules for elements of matrices X and 


is satisfied. Similar to previous considerations, these eigen¬ 
vectors can be arranged into A 3 — Ai conjugated pairs: if | a) is 
an eigenvector, | b) = gr^\a*) is also an orthogonal eigenvector. 
Note that rj is a unitary matrix. 

To derive the effective Hamiltonian, we project Eq. £[]) onto 
the subspace of these eigenvectors. The consideration can be 
done for arbitrary dimension, but for the sake of comprehensi¬ 
bility, we concentrate on this situation with a single channel in 
each lead. In this case, A 3 = 3 and Ai = 1 and we project on 
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line gives rise to a linear splitting of the eigenvalues at fur¬ 
ther increase 6. The eigenvalues approach a linear asymptotics 
given by Eq. d26l) . We illustrate the spectrum in Fig.[9jb). 


V. CONCLUSIONS AND DISCUSSION 


FIG. 9: (Color online) Fine structure of GET in the vicinity of a 
symmetry line with strong SO interaction. S = \h ab \ = \h s j,\ = 

+ \h\ 2 is linear to a phases deviation from the line, (a) Neglect¬ 
ing energy-dependence, (b) The effect of energy-dependence. We set 
JD^D~ = 0.3 and ^ D 2 0 + \Dp = 0.1. 


four orthogonal states \a), \b), \a), \b). For a small phase devia¬ 
tion from the symmetry line, S(tp) « S (^ 0) ) + 6S. The matrix 
element ( 6S ) a p can be represented as 


(0,<a|) T <55 



= i(a\r^rjd + - 2r\5r}\fi) = i2X a p. 


( 24 ) 

Here, 6 is the diagonal matrix of the phase deviations. Note 
that we use {a\t\jti 2 ,\P) = 0 to arrive at this. A matrix repre¬ 
sentation of this correction term X a p for the four bases gives a 
Hamiltonian in the vicinity of the symmetry line 


H - X - 


& • h a b 
ih[) + & ■ li 


- ih {) + & ■ h 

& ' ha), 


(25) 


Here, we use r ! = -gr 2 g. The parameters in block diagonal 
components are defined as h a bp + ih a b ,2 = X a b = 2{a\6\b), 
h a bp = X aa = -Xbb = (a|b|a) - (b\6\b), and the same way for 
h-jy Those in off-diagonal components are h 2 + ih 0 = X„ a = 
-X* bb = < a\6\a ) - (b\6\b), and hi + ih 2 = X~ ba = X* b = (b\5\a). 
The eigenvalues of H are 


± y]hl + L' ± V(/*+ ■ /i-) 2 /4 + (h ■ h + ) 2 + I h 0 h + + hxh~ 

(26) 

with L’ = ( \h \ 2 + | h ab \ 2 + \h ah \ 2 )/2 and h ± = h ab ± h- ah . This is 
similar to Eq. ( IT9l >. however, in the case under consideration, 
all elements of the effective Hamiltonian are proportional to 
the phase deviations. This results in a linear splitting of four- 
hold degenerate eigenvalue and two bands touching the gap 
edge at the symmetry line (Fig.|9ja)). 

The energy-dependence is taken into account the same way 
as Sec. IIV Cl It gives a correction term in the Hamiltonian that 
importantly does not vanish at 6 —> 0, 


// = A + 


D ab D 0 -i& ■ D 


, Do + i& ■ D 


D 


ab 


(27) 


with D a b = (a\D\a) + ( b\D\b) and the same for D ab , Do + iD 3 = 
(a\D\a) + (b\D\b ), and ZU - iD\ = (b\D\a) - (b\D\a). Here, 
D a bD llh > D^ + \D\ 2 due to the causality. This guarantees the 
positive eigenvalues at the symmetry lines, 6-0, given by 

( Dab + D ah )l2 ± yJ(D ab - Dab) 2 14 + D l + \ D \ 2 - The eigenval¬ 
ues are doubly degenerate. The deviation from the symmetry 


We have studied the singularities and peculiarities in the 
ABS spectrum of a Josephson junction connected to N su¬ 
perconducting leads. The ABS energies in such junctions de¬ 
pend on A -1 independent superconducting phase differences, 
E(ipi,<f 2 , ■ ■ ■ ,)fN- 1 ) being a periodic function of all phases. 
Therefore, they can be regarded as energy bands in the N - 1 
dimensional periodic solid, if one associates ipj with quasimo¬ 
menta. We have concentrated on the singularities related to 
topological properties and use numerical illustrations of the 
spectrum and as well as derive effective Hamiltonians to de¬ 
scribe the vicinity of the singularities. The illustrations are 
made for a four-terminal short junction. In this case, the ener¬ 
gies of ABS correspond to the bands in a 3D solid. The ABS 
energies are calculated from Beenakker’s determinant equa¬ 
tion using scattering matrix. 

We reveal the singularities in the vicinity of zero energy 
and near the gap edge. We establish a mathematical analogy 
between the spectrum at E = 0 and E — A. 

First, we have considered Weyl singularities near zero en¬ 
ergy. When the SO interaction is absent, the singularities are 
found at E = 0 accompanying conical spectrum. The Weyl 
singularities occur at isolated 0D points in the 3D space of the 
phases. The SO interaction splits the singular points to mir¬ 
ror symmetric positive and negative energies in spin. A small 
modification of scattering matrix only shifts the position of the 
Weyl singular points but does not eliminate those since they 
are topologically protected. In the presence of SO interaction, 
zero energy points in the vicinity of Weyl singularities form a 
2D manifold in the 3D space of the phases. This 2D manifold 
encircles the singular point. To prove our numerical results, 
we have derived an effective Hamiltonian that is valid in the 
vicinity of the singularity and at weak SO interaction. Eigen¬ 
values of the Hamiltonian reproduce a conical singularity in 
the spectrum. They reach zero at an ellipsoid enclosing the 
singular point. 

Exploiting the mentioned analogy between the ABS at E - 
0 and E - A, we have investigated the spectrum in the vicinity 
of the gap edge. The analogy implies the GET occurs at a 2D 
surface in the 3D space in the absence of SO interaction and at 
isolated points in the presence of SO interaction. Thus, the SO 
interaction generally lifts the GET except specific situations. 
We have established effective Hamiltonians for two specific 
cases: symmetry lines and isolated points. We have also taken 
into consideration a weak energy dependence of the scatter¬ 
ing matrix relevant for fine structure of GET. The effective 
Hamiltonians derived prove the generality of our numerical 
results. For the GET in a strong SO interaction, the Hamilto¬ 
nian indicates a Weyl singularity with a conical spectrum. The 
singularity is shifted from the edge by the energy-dependent 
term similar to the effect of SO interaction at £ » 0. How¬ 
ever, we do not find the isolated GET point in our numerical 







































10 


simulations. 

At a symmetry line, since three of the four phases are equal, 
the four-terminal junction can be regarded as a two-terminal 
junction with unequal channel numbers in the leads. Then, 
only a single ABS is sensitive to the phase differences while 
the other stick to the gap edge. We have also derived an ef¬ 
fective Hamiltonian in the vicinity of the symmetry lines. The 
energy-dependence of the scattering matrix lifts the GET and 
adds another Andreev state, which is localized in the vicinity 
of the symmetry line. 

Our study provides a numerical and analytical evidence for 
the Weyl points in multi-terminal Josephson junctions. These 
points are different from Weyl singular points from that pre¬ 
dicted and recently found in 3D solids, such as TaAs^. - — , 
TaB24, and NbAs^. In the 3D solids, the absence of inver¬ 
sion symmetry in the material and removing the spin degen¬ 
eracy are essential for a stable Weyl point. However, multi¬ 
terminal Josephson junction shows stable Weyl points even 
though the ABS energies are not spin-split. Our study thus fa¬ 
cilitates an alternative method of realization of Weyl fermions 
in condensed matter. 

The Weyl points in multi-terminal Josephson junction 
should provide positive and negative topological charges in 
the space of the phases. The charges are sources of the Berry 
curvature fields. Riwar et ai— have proposed an experimental 
setup to detect the Chern number by measuring the quantized 
transconductance, in similarity with the quantum Hall effect. 
In their scheme, the Chern number is defined as the integral 
of the Berry curvature field on a 2D plane in the space. Phys¬ 
ically, the integration is achieved by sweeping the phases by 
applying finite bias voltages to the superconducting leads. Ri¬ 


war et al. considered a situation without SO interaction where 
the Weyl singular points appear always at E=0 and are dou¬ 
bly degenerate with respect to spin. Then the lowest positive 
(and the highest negative) levels of the ABS are relevant to 
the Chern number. Near the Weyl point, the quantization of 
transconductance may be sensitive to the temperature of the 
junction owing to the undesired thermal activation of a quasi¬ 
particle in the lowest subband. In this paper, we take into 
account the SO interaction and consider Weyl singularities of 
two types. Although the SO interaction shifts the conical point 
from E=0, it does not influence the Chern number arising from 
the singularity. Therefore, the topological charge may be still 
detected by measuring the quantized transconductance. The 
same applies to the singularity at the gap edge. The SO shift 
of the singular point from the Fermi energy implies that the 
probability of undesired thermal activation of a quasiparticle 
state near the point of the singularity is small at sufficiently 
low temperatures, this facilitates the observation of the topo¬ 
logical effect. 
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